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Abstract. In this paper we investigate a nonlinear stochastic partial differ- 
ential equation (spde in short) perturbed by a space-correlated Gaussian noise 
in arbitrary dimension d > 1, with a non-Lipschitz coefficient noisy term. The 
equation studied coincides in one dimension with the stochastic Burgers equa- 
tion. Existence of a weak solution is established through an approximation 
procedure. 



1. Introduction 

Let d > 1 be an integer and D : — [0, l] d . In this paper, we study the following 
nonlinear stochastic partial differential equation on R + x D 

du d ( d . I 2 / \ 

-dt {t ' x) = dz-{dx- U{t ' X)+ 2 U{t ' X) 

+ a(u(t,x))F(t,x), t>0,x€D,i = l,...,d (f.l) 

with Dirichlet boundary conditions 

u(t,x)=0, t>0, xedD, (1.2) 

and the initial condition 

u(0,x) = Uq(x), x € D, (1-3) 

where uo is a continuous function on D with values in the interval [0, 1] and F is a 
noise on R + x D that is white in time and colored in space. The noisy coefficient 
term a is a real-valued function which will be described precisely later. We will 
refer to the problem (fTTTj) - (fOl) by Eq{d, u ,a). 

It is worth noting that the equation Eq(d, ito, 0) gives the classic Burgers equation, 
which arises in physics and have extensively been studied in the literature, see e.g. 
[H] and [23]. 



A class of equations of type Eq(d, uq, c), with noise depending only on time, was 
studied by Gyongy and Rovira [17] . The authors took the space variable in a 
bounded convex domain of M. d , and the class investigated contains a version of the 
one-dimensional Burgers equation [21 HI E [HI [H3 HH] as a special case. They 
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proved existence, uniqueness and a comparison theorem under assumptions on co- 
efficients which include global Lipschitz condition on the diffusion coefficient a. 

However, the present work is concerned with a non-Lipschitz diffusion coefficient. 
Since the F is a space-time noise and a is non-Lipschitz the results of |17| cannot 
be applied here. 

The equation Eq(l,uo,<j) was studied by Kolkovska [2D]- The author proved the 
existence of a weak solution by using an argument of tightness and solving a mar- 
tingal problem. The case considered corresponds to the stepping stone model which 
describes migrating populations consisting of two types when the total population 
size does not change over time. For more details on the genealogy of a variation 
of this model, one can see Chapter 7 in [29J . Our study extends this of [20] to a 
multidimensional case, and with a more general class of diffusion coefficients. 
It is also worth mentioning that a one dimensional version of , perturbed with 
a space-time white noise, was investigated by Adler and Bonnet (TJ [4] on the whole 
real line and with the non-Lipschitz noisy coefficient tr(r) = \pr. The authors es- 
tablished the existence of a weak solution and its Holder continuity. Moreover, they 
discussed (but did not prove) the uniqueness of the solution. 

In the present paper, we deal with a d-dimensional spdes. As it is well known, we 
can no longer use the space-time white noise for the perturbation; indeed, in di- 
mension d > 2, the fundamental solution associated with the operator — A is not 
square integrable when d > 2, while the one associated with the operator -Jpp — A 
becomes less smooth as the dimension increases. Hence, the martingale measure 
approach, introduced by Walsh [30J , to investigate multidimensional SPDEs driven 
by a space-time white noise gives solutions only in the space of random distribu- 
tions (see e.g. [2] [Tfl [26]). 

In order to circumvent this difficulty, Dalang and Frangos |7] suggested to replace 
the space time white noise by some Gaussian noise which is white in time and which 
has some space correlation in order to obtain solutions in the space of real- valued 
stochastic processes. See [SJ and [57] and the references therein for more literature 
in the subject. 

The aim of this paper is to study the existence of solution of a non-linear equation 
Eq(d,uo,<j) in arbitrary dimension d > 1, with a non-Lipschitz coefficient noisy 
term. The idea to prove our main result is to adopt an approach going back to 
Funaki [T3] . 

The paper is organized as follows. In the second section we give a formulation 
of the problem. Section 3 is devoted to define a spatial discretization scheme of 
and obtain a system of stochastic differential equations. The existence and 
uniqueness of a unique strong solution for this system. Section 4 is devoted to 
the tightness of these approximating solutions. The existence of a weak solution 
of Eq(d,uo,a) is established in Section 5 by solving a martingale problem. We 
conclude with an Appendix containing some technical results needed along the 
paper. 
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Note that all real positive constants are denoted by c regardless of their values and 
some of the standing parameters are not mentioned. 

2. Framework 

Let (Q, F, P) be a complete probability space. Let F = {F(<p),<p £ V(R d+1 } be 
a mean-zero L 2 (J7, J 7 , P)-valued Gaussian process with the covariance functional 

/■+OC n r 

E(F(<p)F(rl>))= / da dx dyip(a,x)f{x,y)i)(a,y), (2.1) 
where / : R d x R d — > K denotes the correlation kernel of the noise F. 

Several authors discussed the weaker assumptions in such a way that (|2.1| defines 
a covariance functional, see e.g. the references [71 ISl 129] . 

In this paper, we will take / symmetric and belonging to Cb(D x D), the space of 
continuous bounded functions. Following the same approach in [B], the Gaussian 
process F can be extended to a worthy martingale measure M — {M(t, A) := 
F([0,t] x A) : t > 0, A £ Bb(D)} which shall acts as integrator, in the sense of 
Walsh [3U], where Bb(D) denotes the bounded Borel subsets of D. Let Gt be the 
completion of the cr-field generated by the random variables {M(s,A), < s < 
t, A G Bb(D)}. The properties of F ensure that the process M = {M(t, A), t > 
0, A E Bb(D)}, is a martingale with respect to the filtration {Gt : t > 0}. Then one 
can give a rigorous meaning to solution of the formal equation (jl.ip . A stochastic 
process u : £1 X M + x D — ► M., which is jointly measurable and C/t-adapted, is said to 
be a weak solution to the equation Eq(d, uq, a) if there exists a noise F of the form 
(|2.ip such that for each 93 £ C 2 {D) such that (p = 0on 3D, and a.s. for almost all 
t > and x £ D, 



u(t,x)(p(x)dx = / uo(x)ip(x)dx + / / u(s,x)Aip(x)dxds 

D JD JO JD 



\SJ D u2M ^X {x)dxds+ L 1° 



(u(s, x))(p(x)F(ds, dx). 



In order to formulate our main result, we assume the following hypothesis 



(H) The function a is Holder continuous of order i < a < 1 i.e. there exists a 
positive constant c such that for all x, y £ [0, 1] 

\a(x)-a(y)\<c\x-y\ a . 

Moreover, we assume that a is not identidically zero and satisfies <r(0) = <r(l) = 0. 

Theorem 2.1. Assume that uq is a continuous function on D with values in the 
interval [0,1], and that a satisfies (H). Then there exists a weak solution to the 
equation Eq(d, uq, a). 

Remark 2.1. The hypothesis (H) implies the existence of a positive constant c 
such that for all x £ [0, 1] 

a{x) < cmin(j; Q ,(l -x) a ). 

Examples 
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1. : The function u\ : x i — > 'x(l — x), which corresponds to the stepping 
stone model, satisfies (H) with a = \. 

2. : Let 7 1]. The function a 2 defined on the interval [0, 1] by 02(0) = 
and (72 (x) = —x 1 log(x) for x e]0, 1] satisfies (H). It is a concave function 
with values in [0, 1]. Indeed, 02 is twice continuously differentiable on ]0, 1[ 
and <r%(x) = x"*- 2 {l - 2 7 + 7(1 - 7) log(as)), for all x G]0, 1[. Moreover, 
taking into account the continuity of 02 on [0, 1] we deduce its concavity 
on [0, 1]. It follows that |cr 2 (a;) — (J 2 (y)\ < sup{|CT2(^ — y) — cr 2 (0) | , |<t 2 (1) — 
0-2(1 — x + y)\}, for all x,y £ [0, 1] such that x > y. Finally, we point out 
that \<r 2 (x - y) - cr 2 (0)| = |cr 2 (x - y)\ < c\x - y| 7_e , where e G]0,7 - §[, 
and |<t 2 (1) - cr 2 (l - x + y)\ = |<r 2 (l - as + y)| < c|as - y\ n for any 77 e]0, 1[. 
The result follows. 



3. The approximating processes 

Let n > 2 be an integer. For Ik := (ki,...,kd) £ {1, ...,n — 1} , set x^ = 
^, and define the interval 

J n n 

Consider the discretized version of Eq{d, uq, a): 

1 d 

du n (t, xl) = A n u n (t, xl)dt + - V?(u n (t, xl) 2 )dt + n d a{u n (t, xl))F[dt, H£), 

i=i 

u n (0,xl) =u (xl), 

(3.1) 

where, the operator in the first term is given by 

d 

A n u n (t,xl) = n 2 Y,(u n (t,xl + -) + u n (t,xl~^)-2u n (t,xl)), 
z — ' \ n n I 

i=l 

and in second term by 

VUu n (t,xl) 2 ) = n (u n {t,xl + ^) 2 -u n (t,xl) 2 ), 

with {ei, 1 < i < d} denotes the canonical basis of R d , and 

F(t,K)= I I F{ds,dx). 
Jo Jtli 

The noises {F(t,Jl^),t > 0}, derived from colored noise F on K d with covariancc 
given by (|2.ip . are one dimensional correlated Brownian motions. Moreover, the 
boundedness of / implies that (F(-, 1^))* < ctn~ 2d , where c is a constant depending 
on /. Indeed 

(F(;El)) t = f f f l K (x)f{x,y)l K {y)dsdxdy 

JO JR d JR d 

f(x,y)dxdy 



i ctn- 2d , 
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where c = max{/(x, y); (x, y) G D x D}. 

Now, we proceed as in [21] and we consider the bijection between the grid D„ := 

K = (v> ^. v) . *«0 e {!. •••> which assi S ns t0 each (tT> ■ - v) 

the integer fci + (fc 2 - l)(n- 1) + ... + (fed - l)(n - 

Set JV := (n — and for each i G {1, ..,JV}, let i n be the unique element 
(£, such that i = k x + (k 2 ~ l)(n - 1) + ..."+ (fed - l)(n - l)*" 1 and 

1 < kj < n — 1 for each j G { 1 , . . , d} . 

Finally, let u"(t) denotes u n (t,i n ). Hence, the system (|3.1[) can be written, for 
i = l,...,d, 



N 1 

d<(*) = + ^b h3 u](t) 2 (t)}dt + n d a(u-(t))F(dt,^) 



3 = 1 

<(o) = u (n- 

Let and B^be the N x N matrices such that 

Bn ^ : = (&ij) 1 < i)J -< i v 

For d = 1, it is well known that we have 



(3.2) 



( a i,j)l 



<i,j<N 



and 



and, 



■il = i, 

-2n if j = i, 
otherwise, 



(3.3) 



n if j = i + 1, 
-n if j = i, 
otherwise. 

(d) 



(3.4) 



For d > 2, it is known from [M] that A n can be obtained by induction using a 
relation of recurrence and the fact that A n ^ is known. For the coefficients bij, let 
us first introduce some notations. For an integer I > 1, let 7g be the £ x £ identity 
matrix and set M := (n — l)" - 1 . One can easily show that 







B (d) 



nl 



B 



M 
(d-1) 







\ 



nI M 
B n d ~^ 



nl 



N ■ 



Note that B„ is the (n — 1) x (n — l)-matrix defined by p. 4 
The main result of this section is the following 
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Proposition 3.1. For any integer n > 1 and for any initial random condition 
u n (0) = € [0,1]* fte system 



W 1 

d«?(t) = (*) + gfeiju^t) 2 ]* + n d a«(i))F(^,n^), 

<(o) =er, 

admits a unique strong solution u n (.) £ C([0, +oo[, [0, 1]*). 
Proof. 

We consider first the following system 

N r 



(3.5) 



d«?(t) = £ 

l<(o) = C, 



i. 



dt + h(u?(t))F(dt,IV>), 



(3.6) 



where ft : R — ► R is defined by h(x) — n d o-(x)l{ < x <i}, and g : R — > R is defined 
by g(x) = £ 2 l{-i< x <i}- In order to show that the continuous coefficients of the 
finite dimensional system (|3.6[) satisfy the growth linear condition, we write it under 
the following form 

dU n (t) = A^U n (t)dt + l -B^(U n (t))dt + K(U n (t))dF t 
U n (0)=u n (0), 

with U n (t) := «(t), I < * < N) , S(C/"(t)) := (ff«(t)), l<i<N), K(U n {t)) := 
diag(h(uf(t)),l <i<N), and F t := (F(t, I™„), 1 < % < nY 

Since the first term in the right hand side of (|3.7[) is linear, then it satisfies the 
linear growth condition. For the second one, using an elementary property of matrix 
norms and the equivalence of norms on M. N , we can write 

\\Bi d ^(U n (t))\\ < ||4 d )||||L(C/"(t))|| 

N 

< cJ2\<(t)\<o\\U n (t)\\, 



where the constant c depends on n and d. Concerning the noise coefficient, owing 
to Remark l2.II we obtain a.s. for t > 0, 



N 



N 



5>(«?(t)) a < cn 2d J2<(t) 2a l{0<u m <i} 



i=i 



i=i 

JV 



< en M £(l + «(t)) 2 



i=l 



< cn 2d N(l + \\U n {t)\\ 2 ) . 

Thus, the existence of a weak solution with continuous sample paths for (|3.6[) is 
given by Theorem 16.11 in the Appendix. 
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Now, we will show that < uf(t) < 1 for all i = l,...,n and t > 0. Taking 
p{x) — x 13 with 1 < (3 < 2a, and using the hypothesis (H) and the continuity of 
the trajectories of it™, we get fot < t < T 

Jo 7KW dMs = n I — — ds 

< n 2d T 

< +OG. 

Then, by Lemma l6Tl we obtain that the local time, L^(uf), of the semi-martingale 
m™ satisfies L° (a;™) = for all t > 0, a.s, and i = 1, ■ ■ ■ ,N. 

Moreover, set x~ : = max(0, —x) and by applying the Tanaka formula and summing 
over all indices i = 1, .., N, we get for t > 0, 

N „t N N 1 t N N 

= ^l {K(s)m ^2a l . J u](s)ds~ - / ^l {nr(s) < 0} ^6 ij5 (^ l (s))ds 

i=l Jo i=l j=l ZJo »=i j=i 

:= J1+/2. 

Since — 2; < x~ and taking into account the fact that a^j > for all j ^ i and 

„ t A iV 



h = - ^^(sj^a^l^n^o}^ 

t A iV 

= -/ EK"( s )(" 2 ™ 2l KW<o}+ H ^jJ-WWIO})] 

t A A 
t A A 

= / K'( s )~ (- 2 " 2 ) 1 {« J "M<o} - w"(s) ^ aijl {Kis) < 0} ] 

Jo 3 = 1 i=1M3 
f t N / N \ ft N 

^ / $")«"(*)" SoiJ^W^o} <*s < 2dn 2 / 5]u?(»)-<fa. (3.8) 
J o 3=1 \i=l ) J ° i=l 



For ^2, notice first that by using the definition of g, 

,t A / A 



/a = 



2^ 5I 1 {-i<«?(s)<o} |nd<(s) 2 - h,jUj{s) 2 \ ds. 



=1 \ j=hi¥4 

Taking into account the fact that bij > for all j 7^ i, we infer that 

rt at 



I f 

h < 7; nd / Y] l{-i<<( s )<o}<(s) 2 ds. 



Moreover, for < s < t such that —1 < uf(s) < 0, we have 
<( S ) 2 = (-<( S )) 2 <-<( S ) = <( S )-- 



cS 
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Hence, 



From the estimates 



h<\nd f 2 
jHJ) and (|3.9p , we deduce 

i=i - 70 i=i 



(3.9) 



Then, by Gronwall's Lemma, we obtain that 53<=i u ?(*)~ = 0- Consequently, the 
solution is non- negative for each t > 0. 

Concerning the fact that < 1, it can be easily checked by using the same 

arguments as above for (1 — uf(t)). 

Therefore, the system Q3.5[) admits a weak solution with trajectories in C([0, +oo[, [0, 1] 

To show the pathwise uniqueness for (|3.5[) . we assume that u 1,n (t) : = (u]' n (t), 1 < 
% < N) and u 2 ' n {t) : — (u 2 '"(£), 1 < i < N) are two weak solutions to (|3.5p . with 
the same noise F and the same initial data. 



Set y" 



u) ,n — it 2 '™, then for t > 0, 



N r t 1 AT t 

5>„-/ y^s)ds+-Y J b itj [u)^{s] 2 -u 2 : n {sf)ds 



«7(«J' n (5))- ff (ttJ' n ( fl )) F(d«,]I?.). 



Consequently, by the boundedness of the correlation kernel /, 



* - n 2 
l,n/„\\ _/..2,n 



ds. 



Owing to Remark 12.11 and Lemma 16. 1[ we obtain that a.s. h\ (y™) — 0, for all 
i G {1, N}. Applying the Tanaka's formula for the continuous semimartingale 
y", it follows that 

\Vi(t)\ = / sgn(y?(s))y?(s)ds 

1 N c l 

+ 2 E b ^ J (^"(*) 2 - u- n {s) 2 ) ds 

+n d f ' sgnWis)) \a{u]> n {s)) - a(uf n (s))] F(ds,Tl? n ). 
Jo L J 

Noting that the sign function is defined by 

, , f +1, x>0: 
sqn(x) = < . _ 
y w 1-1, x < 0. 
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By summing over all i G {1, ■■■,N} and taking the expectation, we get 



N \ IN 



E (E = E ( E a ^ I *9"(V?(«M(*)* 

+ 1e [E hi J*sgn(y?(s))(uf n (s) 2 -uf n (s) 2 ) ds 

:= ,h + J 2 . (3.10) 
For Ji, the boundedness of the coefficients a^- yields 

TV \ 



-t JV / AT \ 

j x = e / E s -^(y"( s M( s ) E a ^ 

i=l Vi=l / 

< cE / y>?(«)|d8, ( 3 - n ) 



the coefficients fr^j and the fact that the solutions are with values in the interval 



where c is a constant depending on n. Concerning J2, using the boundedness of 
the coefficients bij 
[0, 1], we can write 

r t N / N 



' h = l E I ^sgn{y^{s))v^{s){u\' n {s) + n 2 ' n {s))[^b i Ads 

Jo j=l Vi=l / 

< cE / ElJ/"(*)l^- (3-12) 

Combining the estimates (|3. 101) — (|3. 12|) and applying Gronwall's Lemma, we obtain 
for all i > 

AT 

EEly?(*)l=o. 

z=l 

Finally, the well-known theorem of Yamada and Watanabe [3T] (see also Corollary 
3.23 in [TH]) implies that (|3.5p has a unique strong solution. □ 



4. The tightness of the approximating processes 

As in [28J, let Y n be a simple random walk whose generator is the discrete Lapla- 
cian A™ on the lattice D^. To get the system (|3.2I) in its mild form we define the 
fundamental solution associated with A™ . 

Fori", j™ G D*, set 

p2(t,i n ,j n )=n d V(Y t n =j n \Y? = i n ). 
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The system Q3.2p can be now written in its variation of constant form 

N , ft N 



(*) - ^E^,i n ,i>o(r) + ^/ E^(* 

i=i J ° i=l 



s ) i n ,j n )b ilj uUsfds 



3 



+ / X)p3(*-*.i n ,rM«?(*)W*,i?»). 

7 ° ,-=1 

In the remainder of this paper, we denote uf (t) by u n (t,i n ), where the relation- 
ship between i and i™ is described in the previous section. Moreover, in order 
to define u n (t, •) for all x £ D, we set K ra (a:) = i n for x € Z™, and p^(t,x,y) — 
Pd(Pi n n {y)) for x,y £ D. Using a polygonal interpolation like as in |28| , we 

can write for all x G D 

u n (t,x) = [ Pa{t,x,K n (y))u {K n (y))dy +1- [ [ p n d (t - s,x,K n (y))u n (s,y) 2 dyds 

JD 1 JO JD 

] Pd(t - s i K n (y))a(u n (s, y))F(ds, dy) 

JO JD 

= v2(t,x) + t$(t,x) + v%(t,x). (4.1) 
The main result of this section is as follows. 

Proposition 4.1. For each T > 0, the sequence {u n ,n G N} is tight in the space 
C([0,T],H), where H: = C (£>, [0, 1]). 

Proof. Note firstly that u n (t, x) are in [0, 1]. Consequently, there exists a constant 
c(p, T) > such that for any p > 1, 



supE 



sup sup \u n (t, x)\ p 

0<t<TxGD 



<c(p,T). (4.2) 



To handle the stochastic integral, we will use the factorization method introduced 
by Da Prato el al. in [10]. Indeed, it is well known that one can write for j3 G (0, 1), 

Pd(t ~ s, x, y)o-{u n (s, y))F(ds, dy) 

(t - sf- x fi(t - s, x, y)Y(s, y)dyds, (4.3) 



JD 

t 



sm((3n 
where 



o JD 



Y(s,y) : = / / (s-r)-f j p n d ( S -r,y,z)o-(u n (r,z))F(dr,dz). 

J0 JD 

In the sequel, we will take < j3 < A. Hence, by Burkholder's inequality, the 
boundedness of a and / and the estimate i) of Lemma I1T21 we obtain for any T > 
and p > 1, that there exists a constant c(T,p) > such that 

sup sup E\Y (s, y)\ p < c(p,T). (4.4) 

«6[0,T] y£D 

Let us now prove an estimate on the increments in both the space and time variables. 
Firstly, we are concerned with the third term in (|4.ip . Clearly, for all x,y £ D and 

< s < t < T 

E|i£(i, x) - u n 3 (s, y)\P<c (EK(i, x) - u£(i, y)\* + E|uJ(t, y) - v%{s, y)\*) . 
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By Burkholder's inequality and the fact that a and the kernel noise / are bounded, 
we can write 



E\u™(t,x)-u%(t,y)\P = E 



JD 



< c 



[Pd(t - s, x, z) - p%(t - s, y, z)] cr(u n (s, z))F(dz, ds) 

2 \ p/2 



(^J d \Pd(t ~ s > x , z ) -Pd(t- s > Vi z )\ dz ^j ds j 



< c\\x — y\\(« 



(4.5) 



where 1 < q < 2. Note that in the last inequality, we have used the estimate it) of 
Lemma 16.21 On the other hand, we have 



W{t,y)-vS(8,v)\* 
ft 



= E 
< c I E 



JD 



[(t - rf-'m - r, y, z)-(s- rf^p^s - r, y, z)l [0>s] (r)] Y(r, z)dzdr 
(s - rf-^p^it - r, y, z) - p%{s - r, y, z))Y(r, z)dzdr 



JD 



E 



E 



{{t - rf- 1 - (s - r y- L )p2(t - r, y, z)Y(r, z)dzdr 



J D 

•t 



(t - rf~ l pl{t - r, y, z)Y(r, z)dzdr 



Is JD 

A t +A 2 + A 3 . 



By Holder's inequality, (I4.4[) and i) of the Lemma RT21 we get 

A 3 < c\t-s\^- xd)p , 



(4.6) 



where < j3 < \ and < A < 4. Concerning A%, the same arguments as above 
yield 



A 2 < c[ I [{s-rf- 1 -{t-ry~ L ]{t-r)- Xd dr 



< c|t-sr 2Adp 

< c\t- sf^-^l 

For ^4i, we use firstly Burkholder's inequality 



[(s-r)^ 1 - {t-rf-^d 



/3-ii 



(4.7) 



Ai = E 
< cE 



(Pd(t - r, y, z) - p%(s - r, y, z))a(u n (r, z))F(dr, dz) 



JD 



(Pd(t ~ r, y, z) - p2(s - r, y, z))a(u n (r, z))f(z, x) 



JD JD 



p/2 



x (Pd(t ~ r, y, x) - p^s - r, y, x))a(u n (r, x))drdzdx 
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and since a and / are bounded, and by iv) of the Lemma 1672 



Ai < c / / / (p2(t-r,y,z)-p n d (s-r,y,z))\z-x\- a (4.8) 
X (Pd(t ~ r,y,x) - j%(s - r,y,x))) p/2 

p/2 

V <{ I ||pg(t -r.//.-l-y/i(.s -r.//. 

< c|t-s| (1 -^)f. (4.9) 

Therefore, taking into account (|4.5l) - (|4.8p . there exists a positive constant c and 
< <5i < | and < S 2 < | such that 

EK(t, y) - y)\ 2p <c(\t- s\ 5 ^ + \\x- y\\ 5 ^) . (4.10) 

From the Proposition 3.2 in |24| . we deduce that (|4.10p holds for m™. Concerning 
«2 , using Holder's inequality, (|4.2H and the estimates on given by the Lemma 
16.21 we deduce that (|4.10l) holds for u^. Finally, we deduce the tightness of the 
sequence {u n ,n > 1} by the Totoki-kolmogorov criterion, see Theorem 4.10 and 
Problem 4.11 in [191. □ 



5. The weak solution 

Since the sequence {u™, n > 1}, is tight in C([0,T],"H), it is relatively compact 

in this space. Thus, there exists a subsequence, that we still denote by u n , which 

converges weakly in C([0,T],'H) to a process u. By Skorohod's Representation 

Theorem, there exists a probability space f2, and on it a sequence u n , as well as a 

v 

noise F and a process u such that («", u, F) = (u n , u, F) and u n converges almost 
surely to u in C ([0,T],H). We will show, by solving the corresponding martingale 
problem, that u is a weak solution of the equation Eq(d,uo,cr). 

Proposition 5.1. For any (p £ C 2 (D), such that if = on 3D, we have 

A4 v (t) : = / u{t, x)ip{x)dx — I uo(x)Lp(x)dx — / / u(s,x)A<p(x)dxds 
Jd J d Jo Jd 



is a martingale with the quadratic variation 

(M<fi}t= / / o-(u(s,x))(p(x)f(x,y)a(u(s,y))(p(y)dxdyds. (5.2) 
Jo Jd Jd 
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Proof. We consider the scheme (|3.1|> and by multiplying both its sides by ^d<p(Xfc) 
and summing over all the elements of the grid D^, we obtain 

71 n 71 n JO 71 n 

X ls X ]k X fc 

= ^ E fin (*> *£M^) - 4 E ^m^e) - f — d E ^ <)A>(^)rf S 

x :k ^ic 

+1 f^EE( 5,, ( s '^)) 2 v,>(4)& 

2 J ° " x£ i=l 



E^kM^M)W^)- 

Note that {Jvt™(t),t > 0} is a martingale as it is a finite sum of martingales. 
Moreover, by (|2.ip and the boundedness of a, we can write 

E(M^(t)) 2 = E / ds I dx I dycp(x)a(u(s,x)My)a(u( S) y))f(x,y) 

„n „n JO Jl? Jl", 



< t V / cfc / dytp(x)f(x,y)ip(y) 



< t dx dy<p(x)f(x,y)<p(y). 

J D J D 



Hence, there exists a constant c depending only on y>, k and T such that 



sup supE(.M™(*)) 2 < c < oo. 



0<t<T n 



Consequently, {A4",n > 1} converges a.s., as n — > oo, to the martingale A4 V 
given by (|5.ip . Also, we have 



71 >OC 



where 

Kit = ( E / / 







E *.MWP(^)Whi)«.(.)Wf^))«* 



which converges a.s. to (|5.2[) as n — > oo. □ 

Proof of the Theorem 12.11 Following Walsh [3U], there exists a martingale mea- 
sure M with quadratic variation 

v(dxdyds) = <r(u(s,x))f(x,y)a(u(s,y))dxdyds, 
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which corresponds to the quadratic variation (M. lp )t- Furthermore, we consider a 
measure martingale N independent of M and characterized by (|2.1[) . and we set 



F{t,<fi) = J J -^^-^ 1 Ms,x)^{o,i}} l p(x)M(dxds) + J l {ii ( s ^ )e{0A}} ip(x)N(dxds). 
Therefore, F(t, (f) satisfies 

(F(',tp)) t = / / tp(x)f(x,y)tp(y)dxdyds. 

JO JDJD 

Hence, because of the form of its quadratic variation, it has the same distribution 
as N, and by Proposition 2.5.7 in [53], we can write 

M v (t) = / / a(u(s,x))ip(x)F(dx,ds), 

JO J D 

therefore, u is a weak solution of the equation Eq(d, uq,<j). □ 

6. Appendix 
Firstly, we recall the Lemma 1.0 in [22J, 

Lemma 6.1. Let Z = {Z(t),t > 0} be a real-valued semi-martingale. Suppose that 
there exists a function p : [0, +oo) — > [0, +oo) such that 

dx 



for all e > 0, and 



o Pi x ) 



-d(Z) s < +oo, 



/O P( Z s) 

for allt > a.s. Then, the local time of Z at zero, L®(Z) is identically zero for all 
t > a.s. 



The following lemma gives some useful estimates satisfied by . 

Lemma 6.2. i) There exists a constant c > such that for any t > and A > 0, 

sup sup Hp" (t, x, -)||i < ct~ x exp(— ct) 

n>l xeD 

ii) There exist constants c > 0, 1 < q < 2 such that for any t > and x, y € D , 
sup / (/ \pj(t- s,x,z) -p%(t- s,y,z)\dz) ds < c\\x - y\\ « "5 ; 

n>l Jo \J D J 

this estimate holds with q — 1 when d = 1 . 

For any T > 0, i/tere exists a constant c > smc/i t/iat /or > 

sup sup sup / \\pd(t — s,x, ■) — p2(t + h — s, x, •) || irfs < chh 

n>l x£D te[0,T] JO 

iv) For any T > 0, there exists a constant c > such that for h > 
sup sup sup / \\p^{t - s,x, •) - pj(t - s + h,x, -)\\ads < ch 1 

n>l iED t£[0,T] JO 
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where < a < 2 A d, and \\-\\ a is defined on an appropriate class of functions ip by 
\\<p\\m= I I W{x)\\x - y\- a W{y)\dxdy 



D J D 



Proof. The estimates i), Hi) and iv) were proved in |24) . that are respectively, 
(A. 16) of the Lemma A. 4, (A. 36) and (A. 38) of the Lemma A. 6. Concerning 
ii), using the fact that p2(t,x,z) — Y\f = iPi{t — s,Xi,Zi), where x — (xi,l < 
i < d) and z = (zt, 1 < i < d), the inequality Yl%=i a i ~ Ilf=i — Ei=i \ a i ~ 

b%\ ll^=i+i \ a j\ ilj=i IM f° r au rea, l numbers a\, a^, bj,, the estimate i) and 
the Cauchy-Schwarz inequality, we can write 

/ (Z>'^* _ s > x > z ) ~ Pdit- s iV^ z )\ d ^j ds 

d d 
Y[f{(t - S,Xj,Zj) - ]]_Pi(t - s, yj , 

et I d i-1 -i 

- / ( TT H-Pi ( £ — Olli x / \Pi(t - s,Xi,Zi)-fl{t- s,yi,Zi)\dzi 
Jo \i=ij=i Jo 



\JD 



dz ds 



< c 



< c 



x II ||p?(t-*,fcv)Hi ds 

2 



E / (^- S )- 2A(d - 1) / 



s,^,^j) s,yj,^)| 2 dzjds. 



On the other hand, using notations of Gyongy [13] and Holder's inequality, we get 

(t _ s) -2X(d-i) / ^(t-s.x^z^-p^t-s^z^dzds 
Jo 

n-1 



E / (t-sr^-VeM-Zfit-^ds^ixi)-^^ 

n—1 / „t \ 1/9 

- C ^\J e M-Ht~s)3 2 )) (j 2 \xi-yi\ 2 Al) 



.i 

n-l 



1 



<c|x i -y l |<! 2 . 

Note that we choose A and q such that 1 < q < 2 and < 2X(d — 1)£ < 1. where 
± + 4 = 1. ^ □ 
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We have also used the following result in Section[3j for the proof see Theorem 3.10 
of Chapter 5 in [12]. 

Theorem 6.1. Consider the SDE 

X t =X + f b(s,X s )ds+ f <j(s,X s )dM s , (6.1) 
Jo Jo 

where t > 0, b : R + x R d — > R d and a : R+ x R d — > R d <g> R d , are Borel continu- 
ous functions, M is a d- dimensional martingale in M\ and X is a d-dimensional 
process. Assume further that the R d -valued random variable Xq is To-measurable 
such that VXq 1 = [/,. If there exists a constant c such that 

||&(*,»)||<c(l + H) and \J2 °?j<t,x?j < 0(1 + 11*11), 

for all t > and x € R d , then for all initial condition probability measure /i on R d , 
there exists a weak solution of the SDE id.l]) . 
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